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Abstract 

We explore the phase structure of a holographic toy model of superfluid states in non-relativistic 
conformal field theories. At low background mass density, we find a familiar second-order 
transition to a superfluid phase at finite temperature. Increasing the chemical potential for the 
probe charge density drives this transition strongly first order as the low-temperature superfluid 
phase merges with a thermodynamically disfavored high-temperature condensed phase. At high 
background mass density, the system reenters the normal phase as the temperature is lowered 
further, hinting at a zero-temperature quantum phase transition as the background density is 
varied. Given the unusual thermodynamics of the background black hole, however, it seems 
likely that the true ground state is another configuration altogether. 



1 Introduction 



Non-relativistic superfluids provide a high-precision laboratory in which to probe many-body physics 
in the extreme quantum regime [J]. In an effort to bring the tools of holography [ , , ] to bear 
on these systems, considerable effort has been devoted to studying non-relativistic deformations of 
relativistic examples^ which enjoy z = 2 scaling [7, 8, 9, 1(J, 11]. Unfortunately, such deformations 
generate highly atypical states in the resulting NRCFT whose thermodynamic and other properties 
are tightly constrained by their relativistic births. In particular, they are in general far from the 
superfluid groundstates of the corresponding systems, for which we currently have no description. 

In this paper we examine certain superfluid states in a holographic NRCFT in a probe approx- 
imation. Our strategy is essentially the same as in the AdS case (see e.g. [12, 1.3]): we study 
an Abelian-Higgs theory in the background of a neutral asymptotically-Schrodinger black hole 
[9, 10, 11] in the probe approximation. Several features of the geometry, however, make the result- 
ing analysis qualitatively different. For example, we are now forced to turn on two components of 
the bulk gauge field: A^, dual to a boundary charge current, and ^4^, dual to a boundary Mass^ 
current. By itself, this is not a big deal. What's surprising given intuition from the relativistic 
case is that the boundary value of this second vector component, Mo = A^\q^ weasels its way into 
the dimension of the boundary order parameter as A = 2 it -y/4 + m? + q^M"^. Specifying the 
boundary NRCFT thus requires not just specifying the bulk matter fields and their interactions, 
but also the asymptotic fall-offs of some of the bulk fields. Similar effects arise in the holographic 
renormalization of the theory, which as usual requires introducing counterterms which depend on 
the boundary operator dimensions; here, these counterterms will explicitly depend on the boundary 
values of some bulk fields, too (see e.g. [14] for a discussion of such effects). 

To build a truly NR superfluid, then, we must generate a condensate for a boundary operator 
with non-zero Mass eigenvalue, M 7^ 0. This is the role of the second component of the gauge 
field - in a gauge where the phase of the condensate is constant, the Mass eigenvalue is simply 
M = —qMo. The boundary value of the second component of the gauge field thus controls the 
breaking of the Mass symmetry in the superfluid phase. 

While the background about which we perturb is a 1-parameter deformation of a relativistic 
example, the superfluid state we find is not, and indeed enjoys quite distinct phenomenology from its 
AdS cousins. Fundamentally, the non-relativistic condensate is characterized by one more quantum 
number than in the relativistic case - the Mass eigenvalue, M, of the order parameter - with the 
NR condensate breaking the symmetry generated by the Mass operator, a key signature of a non- 
relativistic superfluid. As we shall see, this leads to a host of interesting effects in the strongly NR 
regime, including the appearance of a thermodynamically unstable high-temperature condensed 
phase which drives the superconducting transition from 2"'^ order to 1** at a multicritical point, 
the persistence of a condensate even in the absence of a chemical potential for the charge density, 

^ Since the non-relativistic conformal group is a subgroup of the relativistic group in one higher dimension, we can 
construct a non-relativistic conformal field theory (NRCFT) by turning on an operator in a relativistic conformal field 
theory (CFT) which breaks the relativistic group to its non-relativistic subgroup. Taking the operator to be marginal 
in the NRCFT [ >, i>] requires it to be irrelevant in the CFT. Holographically, this corresponds to a 1-parameter 
deformation of the geometry which alters the asymptotic geometry from Anti de Sitter (AdS) , whose isometries form 
the relativistic conformal group, to Schrodinger [7, 8], whose isometries fill out the non-relativistic conformal group. 

^In an NRCFT, each primary operator is characterized by not only a Dimension A but also by its Mass, M, where 
Mass is the name of a central extension M in the NR conformal group. In the case of free fermions, M — Mip^ip, 
hence it is often called the "Number" operator - we prefer "Mass" to disambiguate the various meanings of "number" . 
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and reentrance of the normal phase at low temperatures for sufficiently large background density. 

It is tempting to interpret this re-entrance as signaling a zero temperature quantum phase tran- 
sition as the background mass density is tuned. However, the re-entrant normal state is again the 
simple 1-parameter deformation which we do not expect to be the true equilibrium groundstate, so 
we do not expect this probe analysis to be the end of the story. Meanwhile, it remains possible that 
the system is in fact reentrant for all values of the background Mass density as T — )■ 0, where our 
probe approximation becomes unreliable. Resolving these puzzles, however, requires going beyond 
the truncated probe approximation discussed in this paper; we leave them to future study. 

The plan of the paper is as follows. In Section 2 we quickly describe the basic strategy and 
computational setup, with various details elaborated in Appendices. In Section 3 we explore the 
phenomenology and phase stucture of holographic superfluids outside a Schrodinger black hole (an 
analogous study in the background of a Schrodinger soliton [15] is performed in Appendix A - while 
this is not in the same ensemble as the black hole, it provides an alternate example with surprising 
physics of its own). We close in Section 5 with a summary and list of next steps. 



2 The Setup 



Our basic strategy involves studying an Abelian Higgs system, 
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in the probe limit, — )• oo. Here, / = 1 — r^(7rrri)^/^, K = 1 + r'^Q^ and the metric is given in 
string frame. One can think of this as a rather extreme truncation of the charged Schrodinger black 
hole system [16, 17] where we drop the coupling of the vector to the scalar and massive vector of 
the black hole background, or simply as a holographic toy model. The geometry is controlled by 
two physical parameters, the background mass density, 0, and the temperature, T,with the horizon 
located at the radial coordinate rjj = {ttTO,)~^^^ . 

For spatially homogeneous solutions, we can without loss of generality set ^ = and take 
(J) = (^(r) and A = At{r)dt + A^{r)d(^. In Einstein frame, the equations of motion take the form, 
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2.1 Asymptotic Behavior and the Holographic Dictionary 



Near the boundary at r = 0, the vector components behave as, 

At = IJLQ + pqr"^ + ... , A^ = Mo + pur"^ + ■■■ (2.6) 

where the various . . . represent various (possibly non-normahzable) terms whose coefficients are 
entirely fixed by the equations of motion and the values of these integration constants, /ig, pg, 
Mo and pM-""' As usual, pq represents the chemical potential per unit charge, which effectively sets 
the zero of energy in the boundary theory - the gauge-invariant bulk quantity that becomes the 
boundary hamiltonian acting on the operator dual to the bulk matter field of charge q is {idt + qAt)\ 
at the boundary, for plane waves e~*'^*, this becomes (uj + q pg). Thus, one insertion of the charged 
operator Oe~*'^* costs 6E = {to + qpg). As usual, pg computes the induced charge density. 

It might be tempting to think of Mq as a chemical potential for the Mass operator, M. However, 
this is not quite right - it is a superselection parameter. Recall that, holographically, 

M = P^la = -i{d^ - iqA^)\d , (2.7) 

i.e. M is the boundary value of the gauge-invariant ^-momentum in the bulk. The mass eigenvalue 
of a boundary operator dual to a bulk field with ^-momentum I and charge q is thus M = (i—qMo), 
where Mq = A^\q. Like a chemical potential, Mq sets a bias for the mass M, shifting it away from 
its ^-momentum, i. But the mass in an NRCFT is not a parameter, it is part of the definition 
of the theory. Thus, once we fix gauge in the bulk, different values of A^\q correspond to distinct 
NRCFTs, not to a fixed theory with different background fields turned on. In particular, as we 
will see momentarily, the dimensions of various boundary operators depend on A^\g, an unfamiliar 
effect, pm computes the Mass density coupled to ^4^. Henceforth we fix gauge in the bulk such 
that ^ = and M = -qA^\Q. 

As for our charged scalar, near the boundary at r — )• it behaves as 

^ ~ 0ir^- + (j)2r^+ + ... , (2.8) 

where 

A± = 2 ± ^4 m2 + g2jvf 2 . (2.9) 

(Note that we will occasionally write Ai and A2 for A_ and A_(_, respectively.) In the window 
1 < A- < 2, both components (j)i^2 are normalizable, so we may interpret either of (f)i^2 as the 
vev {O), with the other representing the source J. These two choices correspond to alternate 
quantizations of the boundary NRCFT [7, 18]. We will focus on the choice (O) oc (pi and J (x (j)2 
for reasons which will become clear in the next section. 

^In particular, the leading term for At runs as —2pM log(r-). While formally the dominant term, it is determined 
by the equations of motion and pM and thus does not represent an independent mode of the system. Importantly, due 
to factors of the inverse metric, this log running does not lead any components of the bulk stress tensor to diverge. 
A complete holographic renormalization of this system would settle the dictionary, but is beyond the scope of the 
present paper; for the moment we simply take the above dictionary as a provisional interpretation which is supported 
by the consistency of the results below. Interestingly, while At has no log in fully backreacted charged-black-hole 
solutions [16, 17], linearizing the Maxwell equation around these solutions does generate a log without changing any 
other of the asymptotics of the vector, so this log is likely a simple consequence of an extreme truncation of the full 
charged black hole system. It would be interesting to study the full system and see what, if anything, changes. 
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Importantly, the dimensions, A±, depend not only on the mass of the bulk scalar, as in AdS, 
but also on the boundary value of a bulk field. Mo = A^^Iq. As discussed above, by the holographic 
relation M = —i{d^ — iqA^)\Q, this quantity is nothing but the Mass eigenvalue of the dual operator, 

M = -qA^le = -q Mo. 

Our expression for the dimensions above then becomes, A± = 2 it \/4 + + , which is the 
expected form [7, 8], including the quadratic dependence on M inside the radical. 

Now, as discussed in [9, 10, 11], the free energy of the full system takes the form, F = E + ^xmM, 
where /j-m = ^4 is determined by the background spacetime. The total free energy per insertion 
of an operator dual to a bulk field with ^-momentum i, frequency uj, and coupled with charge q to 
our gauge field is thus 6F = (w + q fig) + hmM, where M = {i — qMo) 

2.2 Near-Horizon Behavior and Setting Up the Calculation 

In the bulk, we are thus left with a six-parameter family of solutions labeled by sources (/xg, Mq, J) 
and responses {pq, pM, (C))- Holographically, we expect boundary conditions at the horizon, where 
the radial equations of motion degenerate, to impose three additional constraints. Together with 
the two parameters T and O of the background geometry, this should leave us with a five-parameter 
phase space. To verify this, we need to study the behaviour of our solutions near the horizon. 

The equations of motion degenerate at the black hole horizon, so we must impose boundary 
conditions to pick the appropriate solutions. As usual, it suffices to impose regularity at the 
horizon, which is in any case necessary for the validity of the probe approximation. Assuming 
regularity, the equations of motion as presented in (2.3) - (2.5) degenerate into three algebraic 
equations relating the six horizon values of the fields and their derivatives, as expected, 

(A^irn) -2rj,n^At{rH)f ^{rn) = (2.10) 
-2r|,^;(rH)+(^f^)4(rH) + 2(zVH</'c/>'(r^^)^(r^,) = (2.11) 

ArH(p'{rH)+(4KH{KH-l)rU^Aj{rH)+m^K^/^^^{rH) = 0. (2.12) 

This suggests a simple numerical strategy for constructing superfiuid states of our holographic 
NRCFT. To specify a solution to the full equations of motion, we fix any three of fiq, Mo, J , pq, pm 
and (O) at the boundary and impose the above regularity conditions at the horizon. Since we are 
interested in spontaneously generated condensates, we will generally set J = 0. The resulting two- 
point boundary value problem can be solved numerically in various ways. The most straightforward 
is a brute-force shooting method, as typically employed in the relativistic case. 

In sweeping out parameter space, however, we must be careful to vary the parameters of the 
NRCFT while holding the NRCFT itself fixed - i.e., while holding the spectrum of quantum 
numbers fixed. This is straightforward in AdS, where fixing the set of dimensions reduces to fixing 
the bulk mass m? of the bulk scalar. Here, however, the dimension A and Mass M of the boundary 
scalar operator depend on the asymptotic value of A^ as M = —q A^\q and A± = 2±y/4 + -|- M^. 
Before sweeping out parameter space, then, we must fix A^\q=Mo. As we've already set J = 0, 
fixing the system thus leaves us with a three-parameter phase space labeled by /ig, 17 and T. 
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This peculiar behavior - that the definition of the boundary CFT depends on the boundary 
behavior of the bulk fields - is a very general phenomenon in Schrodinger holography. Indeed, 
renormalizing the boundary stress tensor, say, or other operators in the boundary NRCFT, requires 
counterterms which are local in time and space, but which depend explicitly on dimensions, A, 
and thus on the asymptotic values of A^, in a mildly non-local fashion. Such field-dependent 
counterterms have appeared previously in attempts to renormalize holographic NRCFTs, most 
recently in [14]. A complete understanding of the holographic renormalization of these theories is 
clearly of considerable interest. 



2.3 Conductivity 

We can compute the conductivity in our superconducting background by studying linear response 
to a time-dependent vector potential Ax ~ e~*'^*. As usual, this boils down to solving the equation 
of motion for the bulk gauge component Ax linearized about the superfluid background and subject 
to infalling boundary conditions at the horizon. Setting Ax = a(x)e~*'^*, we have 

fr^a-i - (4 - /^) ™; + (^^^^ - 2,W) = . ,2.13) 



Near the horizon, this reduces to. 



where r = rn — The infalling solutions thus takes the form, 

Ax = aoe-'''\r-rH)^il + ai{r-rH) + ...) (2.15) 

Near the boundary, 

Ax = Ao + A2— + ... 
A short computation then verifies that the conductivity is given by, 

, ^ (Jx) . (Jx) . A2 



{Ex) u:{Ax) ojAo 

Note that, since we are solving a linear equation but only care about this ratio, the overall scale of 
Ax is immaterial. We can use this freedom to set oq = 1, which simplifies the numerical problem. 

Notably, we can analytically determine the w-dependence of a for large and small oj via standard 
power-series analysis. Importantly, the scaling in the superfluid phase to be independent of A and 
M. At small frequency, we find, 

Im[(j(cj < 1)] oc (2.16) 

while for large oj we have, 

Re[a{uj > 1)] oc uj-^''^ Im[a{oj > 1)] oc oj-^/^ (2.17) 

This last result unsurprisingly differs from the AdS case, where i?e[cr(tJ ^ 1)] = 1. Reassuringly, 
they both match numerical results presented below, a nice sanity check. 
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3 Phases of a Schrodinger Superfluid 



Thus armed, we now get down to the business of finding a superfluid state in our non-relativistic 
holographic CFT and exploring its phase diagram. A priori, the phase space is fairly high- 
dimensional - specifying a point involves fixing A and AI to fix the theory, then tuning /ig, T 
and Q to sweep out the phase diagram. For simplicity, we will begin by picking convenient values 
A = 6/5, M = l/2 and /ig = 1/8, then dial the background Mass density, fi. This will reveal a 
zero-temperature quantum phase transition at a critical value 0^,. We will then fix $7 and vary /ig, 
which will drive the superconducting phase transition from 2"*^ to 1'^* order. 

3.1 Varying and a Quantum Phase Transition? 

We begin by fixing A = 6/5 and M = 1/2, then set /xg = 1/8 and vary $7 between and 1. The 
basic results are presented in Figures 1, 2 and 3, which plot the condensate {0{T)) as a function 
of temperature, as well as the AC conductivities Re[a{uj)\ and Im[a{u})], for J] = j^, | and 1, 
respectively. These results are discussed in detail below. 

For very small Vt, the geometry remains essentially AdS until very close to the boundary, so we 
expect most low-energy physics - such as superfluid condensation - to very closely track familiar 
AdS results. This turns out to be almost correct, modulo a surprise we'll explore shortly. 



{0,> 




Figure 1: At small Q., the behavior of the superfluid is essentially the same as in AdS, with a 2 
order mean-field phase transition at the onset of superconductivity at Tc, including the familiar 
gap-and-pole form in the AC conductivity, leaving us in a happy superfluid state at T = 0. Here, 
(/xg, Q) = (1/8, 1/16), with Tc = 0.505. 

Figure la shows the condensate as a function of temperature for the first conformal family, for J7 = 
1/16. As is clear by eye and can be checked precisely from the numerics, the resulting condensate 
turns on at r = Tc with classic mean- field behavior (/3c = \) and grows as the temperature 
is lowered. Figures lb and Ic then show the real and imaginary parts of the AC conductivity 
for various temperatures indicated by color, from high (violet) to low (red). These demonstrate 
the appearance of a superconducting state at T^ with the gap growing as the temperature is 
lowered. Note, too, that the conductivity in the superfluid phase has Im[(j{io — )• 0)] ~ 1/u;, while 
Re[a{uj — )• (X))] ~ Im[a{uj — )• oo)] ~ w"^/^. This scaling is expected on general grounds, so gives us 
confidence in our numerical results. 
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As we increase the background number density, the story changes dramaticahy. Figure 2 shows 
the same plots as Figure 1 but with = 1 rather than Q = 1/16. The most obvious difference 
is that the order parameter vanishes at sufficiently low temperature, T < Tl, doing so again with 
mean-field behavior. As is clear form the finite value of i2e[fT(0)], the extreme low-temperature 
phase is again metallic. 




Figure 2: At large 0,, in addition to the original transition to a superconducting state at Tc, the 
system now exhibits reentrance of the normal phase at a new low-temperature 2"*^ order transition 
at Tl, again with mean-field exponents. Below Tl, a behaves like the normal gas. Here {fXQ,Q) = 
(1/8, 1) with Tc = 0.149 and Tl = 0.009 

Consider now the behavior of the system at zero temperature as a function of the background 
number density, il.. As O — )• 0, the system is superconducting. As — t- 1, the metallic phase is reen- 
trant. At some critical 0=^, then, the zero-temperature system appears to undergo a superconductor- 
metal quantum phase transition. 

It is tempting to try to determine what happens as we tune towards this critical ri*. Figure 3 
shows the same system at O slightly above ri,,,. As before, there is a phase transition at Tc with 




Figure 3: At intermediate we again have a 2 order mean-field transition into a superfluid 
state at Tc- At low temperatures, however, the system undergoes a non-mean-field transition to an 
apparently insulating state. Here (/iQ,r2) = (1/8,3/8) with Tc = 0.123. 

standard mean-field behavior. The zero temperature behavior, however, differs dramatically from 
mean-field expectations; rather, at low temperature, the condensate decays exponentially, as does 
the superfluid density, while the normal density remains vanishing and the conductivity heavily 
suppressed at small but non-vanishing u, suggesting that the T = state is not metallic. It is 
tempting to read this as indicating a translationally-invariant insulating phase. 
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However, numerical results in this region should be taken with a sizeable grain of salt. Indeed, 
at sufficiently low temperature, the numerics simply fail to converge. More physically, in this 
regime, the probe approximation is becoming dangerously unreliable - the matter field profiles 
which generate the required boundary values grow rapidly deep in the bulk (and in particular 
near the horizon) as we approach T = or ^l^,. Backreaction may thus qualitatively alter the 
low-temperature physics, either near the transition at O ~ 17* or for sufficiently low T at any fi. 

Indeed, it is entirely possible that the backreacted solution is re-entrant at any value of fl; our 
analysis is only reliable sufficiently far away from T = 0. To unambiguously exclude re-entrance at 
small 17 as T — )• requires including backreaction, which is beyond the scope of this paper. Note, 
however, that the probe approximation shows no signs of inconsistency for 17 > 17*, so we can be 
quite confident that the system is definitely re-entrant at sufficiently large 17. 

3.2 High Temperature Condensates and the Free Energy 

The surprise alluded to above involves the /li^/i-temperature limit. Figure 4 shows the same system 
but now extending to higher temperatures. The surprise is the appearance of a high tempera- 
ture condensate at T > Th- Troublingly, the condensate appears to grow without bound as the 
temperature increases.^ 



0.12 
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0.08 
<0i) 0.06 
0.04 
0.02 
0.00 

1 2 3 4 5 
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Figure 4: Surprisingly, there is another condensed phase at high temperatures. Here (/iQ,17) = 
(1/8,1/16). 

Before we panic, however, we should verify that this high-temperature condensate is in fact 
thermodynamically favored over the trivial vacuum. Holographically, this means we computing 
the holographically renormalized on-shell action. Unfortunately, in asymptotically Schrodinger 
spacetimes, holographically renormalizing the action is exceedingly complicated. Happily, a simple 
strategy allows us to compute the difference in free energy between condensed and vacuum states 
without performing a full renormalization of the action:'' 

The basic idea goes as follows. Generally, specifying the non-normalizable (source) mode (pi of 
the bulk scalar determines the normalizable (response) mode, (j)2- Smoothly varying the source thus 

''Such a high-temperature instabiUty was predicted by Cremonesi et al whenever A < 4. 

^We thank Nabil Iqbal for iUuminating discussions on this topic, and refer the reader to [20], in which this approach 
is further developed. 
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traces out a curve (/>2(</'i) in the (01, ^2) plane. Along this flow, we can ask how the free energy 
- aka, the Euclidean action - varies. Given the properly renormalized action, the variation of the 
full bulk action takes the form, SSeff = . . . 5(l)i + . . . 5Ai where 6(j)i and 6Ai are the variations of 
the bulk fields and the . . . correspond to the bulk equations of motion. So long as we satisfy the 
bulk equations of motion, this reduces to a simple boundary term, 5S = (Ai — A2) Jq^^ 4>2 S4>i — 
2 Jq^{pm S^q + pq 5Mo) Moreover, if we hold fixed the asymptotic values of Ai (corresponding to 
fixing the values of the chemical potential fiQ=At\d and the mass M=A^\q), this further simplifies 
to, 5S = (Ai — A2) Jq^ (j)2 6(pi We can thus compute the relative free energy density [Fa—J^b) 
between any two states A and B connected by such a flow by integrating 5S along the flow, 



Fb - Fa = - 




where Vd is the volume of the boundary theory and the integral is performed along the flow specified 
above. By construction, this agrees with what we would get by evaluating the fully holographically 
renormalized free energy for each solution and subtracting. Happily, this allows us to compute the 
correct free energy without having to worry about the full holographic renormalization of the theory 
(for further comparison between holographic renormalization and our method, see [14],[2J],[22]). 

Now consider the case of our holographic superfiuid in alternate quantization, where 4>i={0) is 
the response and 4)2= J is the source. In this case, the curve </>i(02) is multi- valued over </>2=0, 
with one solution corresponding to the trivial vacuum, (O)=0, and one to the nontrivial condensate, 
(0)7^0. As outlined above, these two solutions are connected by a very specific flow in the ((/)i, 42) 
plane. To compute the properly renormalized relative free energy, then, all we must do is find this 
fiow and integrate along it, 

Fc-Fn = -T(Ai - A2) (P2 d<Pi (3.2) 

Jn 

where the integration is again along the flow defined above. If this difference is negative, the 
condensate is thermodynamically favored. 

Figure 5 plots two such flows. On the left we have a flow connecting the trivial vacuum ((/>i=(^2=0) 
and a non-trivial vacuum (^17^0, 02=0) of the first conformal family in the low-temperature regime, 
with the flow indicated by the solid line and the direction of flow deflning the direction of integration. 
The area under the curve, corresponding to the free energy of the condensed state, is negative. On 
the right is the analogous flow in the high-temperature regime - here the free energy is positive. 
We thus deduce that the low-temperature condensate is thermodynamically stable, while the high- 
temperature condensate is unstable, at least for this first conformal family. 

What about the second conformal family? Figure 6a plots the condensate of the second family as 
a function of temperature. Note that there is no separate low vs high temperature condensate, just a 
single continuous instability whose profile grows with temperature. Figure 6b then shows a typical 
fiow at typical temperature. Importantly, the enclosed area is negative for every temperature, 
indicating a thermodynamic instability even at arbitrarily high temperature. This is why we quietly 
chose the first conformal family in Section 2. It would be interesting to understand the meaning of 
this instability in detail. 
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Figure 5: (a) Low T condensate has smaller free energy than non-condensed phase, Tq — < 0. 
(h) High T condensate has larger free energy than non-condensed phase, J-q — J^n > 0. One can 
determine the sign of — from the orientation of the curve. Here {fiQ,Q) = (3/8, 1/16). 
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Figure 6: (a) Condensate as a function of T for the second conformal family with A2 = 13/5. (h) 
Typical flow at generic temperature, indicating a thermodynamic instability at every temperature. 



3.3 Varying /ig and a Multicritical Point 

The thermodynamic instability of the high-temperature condensate leads to an important physical 
effect as we vary ^q. Figure 7a plots Tc{^q) and Th{hq), the critical temperatures for the low- 
and high-temperature condensates as a function of the chemical potential /Uq. As we crank up fit, 
holding all other parameters fixed, Tc increases while Th decreases. At a critical value, /U^,, the 
two critical points merge; above ^=1,, the condensate is non-zero for all temperatures. This is clear 
from Figure 7b, where we plot the order parameter as a function of temperature for values of the 
chemical potential above and below this critical ^* . 

However, we have already checked that the condensed phase is thermodynamically disfavored at 
high temperatures. For jiq > ij,^:, then, there must be a critical temperature, T^:, above which the 
smoothly varying, non-vanishing condensate becomes thermodynamically disfavored. This temper- 
ature is indicated by the red dashed curves in Figures 7a and 7b. 

We can verify this by computing the relative free energy of the condensed phase as we vary the 
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Figure 7: At hq > ~ 0.192, the transition goes 1^^ order. Here, 0=1/16. 
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Figure 8: (a) Flow lines in the neighborhood of the critical temperature, for fit = 3. (b) The flow 
switches direction discontinuously (lower solid to upper dashed curves) at a critical temperature T^, 
indicated by the red dashed curve in Figure 7, leading to a first order phase transition at . Here, 
fj,Q=3/8 and il=l/16. 

temperature. Figure 8a shows the flows associated to points to the left and right of the critical 
temperature where the low- and high-temperature instabilities meet (indicated by the red dashed 
curve in Figure 7b). Figure 8b focuses in on the immediate neighborhood of the transition tem- 
perature for fiQ > fi^. For all temperatures below the critical temperature, the flows go below the 
horizontal axis, corresponding to a negative free energy and a thermodynamically stable conden- 
sate. For all temperatures above the critical temperature, the flows go above the horizontal axis, 
so the condensate is thermodynamically disfavored at high temperatures. Indeed, while the value 
of the condensate is non- vanishing and in fact completely smooth as we flow through T^,, the path 
which carries us from the trivial state to the condensate changes discontinuously as we pass through 
T*. As a result, the integrated area - and thus the free energy - also changes discontinuously at 
r=K. Moreover, as we take /xg — )• the value of the condensate at the transition goes to zero, 
(0(T=k)) — )• 0; this ensures that the latent heat of the transition goes to zero at the multicritical 
point where the transition switches from 1^* to 2""^ order, as expected on general grounds. 

The upshot of all of the above is that as we raise fiq , the phase transition from high-temperature 
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metal to low-temperature superfluid switches from 2 order to 1**, with the transition occurring 
at a multicritical point where the low- and high-temperature superfluid phases collide. Near the 
phase transition boundaries, including the multicritical point, the order parameter scales with 
simple mean-field exponents. More precisely, near the finite temperature 2"''^ order phase transition, 
{O) ~ (rc-r)V2, while near the 1^* order phase transition boundary when /ig > /i* the condensate 
jumps discontinuously at T^,^ with {0{T^)) ~ (/ig — yL^,)^/"^. This can be succinctly encoded in a 
simple mean-field free energy, F{(p) = \c2{T — Tc{fiQ))(p'^ + |c4(;U=k — fiq)'^'^ + gCgV'^, with tp ~ {O) 
and with coefficients C2, C4, ce > 0. 

3.4 Setting yUg = and the Persistence of Condensates 

Playing with raises another interesting point. Fundamental to our construction is that the 
scalar operator carries a charge q under a global symmetry of the boundary theory, /xq tells us the 
energy cost for adding a unit of this charge to the system. In AdS, superfluid condensation is often 
induced by tuning /ig beyond a threshold. Is this also necessary in the non-relativistic case? 




0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.000 0.002 0.004 0.006 0.008 0.010 

^^Q T 



Figure 9: (a) Fixing T ^ 1 and varying /it- (h) Fixing fit = and varying T. Here, VL=1/1Q. 

By varying /xg to zero while holding all other parameters fixed (see Figure 9a), we see that 
condensation persists even at vanishing /xg. Indeed, by plotting {0(T)) for ^g = (see Figure 9b), 
we see that the form of this curve is quite similar to the large- case studied above with ^g 7^ 0, 
modulo an overall scaling of the condensate. It is tempting to speculate that this indicates two 
distinct pairing mechanisms, one involving the charge and one involving the Mass eigenvalue alone. 
It would be interesting to explore this point further. 



4 Conclusions and Open Questions 

In this paper we have constructed a toy model of superfluid states in holographic NRCFTs and 
studied the resulting phase diagram, finding several unanticipated features. First, as we lower the 
temperature in the disordered metallic state, the system generally undergoes a phase transition to 
a superfluid state. At small (and even vanishing) chemical potential, this transition is 2"*^ order 
with mean-field exponents; at large chemical potential, however, the transition runs strongly 1^* 
order. Secondly, for large background mass density, the superfluid state only appears in a finite 
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temperature window, with the metahic state reentering at sufficiently low temperature. Finally, 
at zero temperature, the reentrance of the metallic phase leads to an apparent quantum phase 
transition from superconducting to metallic as the background mass density is varied. 

Several features of our results deserve further scrutiny. First, our low temperature results derive 
from a probe analysis which is not valid at zero temperature ~ indeed, as we push the temperature 
to zero near or below the putative quantum phase transition at 17*, the bulk profiles of various fields 
grow rapidly, diverging as we approach zero temperature. To be sure, we checked the consistency 
of the probe approximation in each calculation presented above. However, it is entirely possible 
that various of our results could change qualitatively when we include backreaction. To nail down 
the T = physics, we must incorporate backreaction. 

Relatedly, we have tacitly assumed that the neutral black hole geometry is the dominant saddle 
at r = 0. However, for a variety of reasons including the strange thermodynamics of this black 
hole, this seems unlikely to be the case. It is tempting to speculate that the low-temperature 
phase is dominated by a Schrodinger soliton analogous to the AdS soliton which dominates the 
relativistic case a la Hawking-Page. Indeed, a simple such Schrodinger soliton solution is known, 
and we repeat the above analysis for this geometry in an Appendix. However, the black hole and 
the soliton enjoy incommensurate asymptotic periodicity conditions, so cannot contribute to the 
same ensembles. Understanding the true low-temperature ground state of the neutral black hole, 
even in the absence of any charge density in the system, is of considerable interest. 

To this end, it is worth emphasizing that the basic trouble with the thermodynamics of this - and 
indeed all known - asymptotically Schrodinger black holes is the light-cone relation between the 
near-horizon killing vector dr and the asymptotic timelike killing vector, dt- In most constructions, 
this follows from the structure of the salient solution generating technique. The challenge, then, is 
to build solutions which do not flow to AdS black holes near the horizon. 

Meanwhile, it's important to keep in mind that the simple Abelian-Higgs theory we study is 
an extremely stripped down toy model for which we do not have an explicit charged black hole 
solution. In the few examples where such a solution is known [16, 17], the matter sectors are 
considerably more complicated, which is why we worked with the toy model at hand as a first step. 
It would be interesting to repeat our analysis in one of these more elaborate systems to disentangle 
the peculiarities of our toy model form general features of non-relativistic holographic superfiuids. 

Finally, several intriguing features of this system still need interpretation. Why does the second 
conformal family have a thermodynamically dominant high-temperature instability, what does that 
instability signal, and is there a simple low-energy way to see that this family is disfavored? Our 
conductivity calculations reveal a number of quasiparticle peaks with rather peculiar behavior, 
particularly near the critical point at Vl^ - are these artifacts of the probe approximation, or do they 
signal real physics, and if so, what are they telling us? Does a proper holographic renormalization 
of the full system (which remains an open problem) alter any of our results, and if so how? What 
about the pairing mechanism - what is the fermion spectral function in these systems, and can 
we correlate pairing of probe fermions with the condensation seen above? We hope to return to 
various of these questions in the future. 
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A Superfluids in a Schrodinger Soliton 

It is entertaining to apply the approach developed above to a slightly different geometry, the so- 
called Schrodinger soliton: [ ]: 

dSsolUon,Str = [-fs + J ^ + + -^d^ + 7^ + ^ ' ^""-'^ 

where Kg = 1 — r^O^ and = 1 — r^/r^, with controlling both the gap and the radius of the 
^-direction, L^=-^. Here, the radial direction is cut off smoothly by a spacelike circle shrinking 
rather than by a black hole horizon - and indeed this solution was obtained by double wick rotating 
the Schrodinger black hole with compact ^-direction (which is spacelike near the horizon).*^ 

As in the black hole case, a probe superfluid in the soliton geometry is characterized by five 
parameters: two define the theory (the dimension A and Mass M), two are properties of the 
background which fix thermodynamics quantities (the mass density and mass gap of the soliton, 
mG=jj^) and finally the ?7(l)-charge chemical potential /ig determined by the non-normalizable 
mode of the bulk gauge field At. In the remainder of this appendix we briefly summarize the results. 

A.l Varying /ig 

Figure 10(a) shows the condensate as a function of the chemical potential /xg, revealing a critical 
minimum value at which the system undergoes a 2""^ order transition with mean-field exponent. 
For /i < ;Uc, Im[a{oj — t- 0)] — t- finite, indicating a translationally invariant insulating phase. For 
> /^c! by contrast, /m[cr(a; — )• 0)] ~ ^, indicating superconductivity. We would thus appear to 
find a 2"'^ order insulator-superconductor quantum phase transition by varying //g. However, in 
addition to this w — )• pole, we find two more mysterious poles at cji and uj2 separated by a finite 
gap. This is reminiscent of the paired poles we found in the black hole system at intermediate 
values of the background density near the critical point at ri*. This transition also recalls the AdS 
transitions studied in [ ]. 

®It is tempting to identify this solution as a low-temperature confined phase of the Schrodinger black hole studied 
above, analogous to the Hawking-Page transition from AdS black hole to AdS-soliton. However, this is not correct: 
regularity of the euclidean solutions imposes incompatable periodicity conditions. More precisely, for the black hole, 
smoothness of the global euclidean geometry and compactness of the direction dual to the Mass operator require the 
periodicities itt = itt + n/T and = £.b + nn/T + wL^. For the soliton, on the other hand, we need its = its —in/Ts 
and ^s = ^s + inn/Ts + wL^. The inequivalence of these conditions (note in particular the extra factor of i in Ts) 
tells us that these solutions correspond to inequivalent ensembles. 
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Figure 10: (a) We find a 2"^^ order phase transition at fic = 1-17. (b) When < ^c, Im[a{uj — t- 0)] 
is finite, indicating an insulating phase, (c) For > fic, we find a superconducting ^ pole, as well 
as two gapped poles at finite oj. 



A. 2 Varying Vt 

In the black hole case, tuning Q. drove us through a superfluid-conductor phase transition at zero 
temperature. The (zero-temperature) soliton shows the same effect, with the spontaneous conden- 
sate disappearing in a 2'^'^-order transition as Vl passes through a critical value, Q.c ~ 0.163, as 
shown in Figure 11a. Here, however, the normal phase is an insulator (cf Figure lib) with the 
gap controlled by {Q. — J7c) and a double-pole structure as seen above. Interestingly, this double- 
pole structure persists into the superconducting phase, with the gapped poles merging with the 
zero- frequency poles at a finite value of ~ \^c- 

0.4 : 

0.3 - \ - 

<0|> \ 

0.2 - \ ^ 

0.1- \ ; 

0.0 [ ' . .- 

0.00 0.05 0.10 0.15 

(a) a (b) 

Figure 11: (a) A 2"*^ order superfluid-insulator quantum phase transition at $7c = 0.163. (b) For 
> we find an insulating gap. For ^^Ic < ^ < ^c, we find a superconducting pole, plus two 
isolated poles. For Q < ^f^c; these poles merge with the zero-frequency superconducting pole. 




A. 3 Varying rric and the Gap 

In the above we have held mo fixed. As it turns out, the only effect of varying mc is to rescale the 
gaps in all of the above (see Figure 12). This fits nicely with naive intuition for the effect of the 
compactification radius of the ^-direction. 
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Figure 12: Varying niQ resettles the insulating gttp. Here, mG=^ (left), (right). 

B Scaling Symmetries 

The system described by (2.1) and (2.2) enjoys three distinct scahng symmetries which we can 
use to fix various parameters to convenient values. Re-introducing Gn and Ra in the action and 
metric, the three scahng symmetries act as: 
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(B.l) 

-ic ^ = 1 and Ra = 1. The third, which is the basic 
scahng symmetry of the Schrodinger system with dynamical exponent z = 2, can be used to fix rjj 
to a convenient reference value, r^. Given that T = ^ , this fixes a relation between T and ft. To 



The first two symmetries can be used to fix 



access more general values of these parameters, corresponding to {r'j^,T' we simply map the 
system to (ro,T, il) = (A~^r^, A^T', A^fi'), with all physical parameters correspondingly rescaled. 
(As usual, taking q — )• oo while holding g$ fixed for each matter field $ gives us a probe limit in 
which backreaction is negligible [24].) 
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